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Abstract 

Calculations of the expected TeV 7-ray emission, produced by accelerated cosmic 
rays (CRs) in nuclear collisions, from supernova remnants evolving in a uniform in- 
terstellar medium (ISM) are presented. The aim is to study the sensitivity of 7-ray 
production to a physical parameter set. Apart from its general proportionality to 
Njj, it is shown that the 7-ray production essentially depends upon the ratio of the 
CR diffusion coefficient n to a critical value k ct u = 10(i?o/5 ^G)(-/V#/0.3 cm -3 ) -1 / 3 ^, 
where Bq and Nh are the magnetic field and hydrogen number density of the ISM, 
and Kb denotes the Bohm diffusion coefficient. If k is of the same order or lower 
than K cr it, then the peak TeV 7-ray flux in the Sedov evolutionary phase, normal- 
ized to a distance of 1 kpc, is about 10~ 10 (./V#/0.3 cm -3 ) photons cm -2 s _1 . For a 
CR diffusion coefficient that is significantly larger than n C rit, the CR cutoff energy 
is less than 10 TeV and the expected 7-ray flux at 1 TeV is considerably below the 
presently detectable level of 10 -12 photons cm" 2 s _1 . The same is of course true for 
a SNR in the rarified, so-called hot ISM. 



1 Introduction 



A direct empirical test, whether or not the observed Galactic cosmic rays 
(CRs) are indeed produced in supernova remnants (SNRs) at least up to an 
energy of about 10 14 eV (see e.g. [1-3] for reviews of CR acceleration theory) 
should be possible with observations of SNRs in high energy (e 7 ~ 1 TeV) 
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7-rays. If the dominant nuclear component of the CRs is strongly enhanced 
inside SNRs, then through hadronic collisions, leading to pion production and 
subsequent decay, 7-rays will be strongly produced. 

Theoretical estimates of the 7r°-decay 7-ray luminosity of SNRs [4-6] have led 
to the conclusion that the expected TeV 7-ray flux from nearby SNRs in high 
enough ambient densities should be just detectable by present instruments. 
The conclusions were based on a hydrodynamic approximation for the CR 
component and involved (very reasonable) assumptions about the CR energy 
spectra. 

The spectra of 7r°-decay 7-rays, produced by shock accelerated CRs in SNRs 
that expand into a the uniform interstellar medium (ISM), were studied in 
detail in a previous paper [7] (hereafter referred to as Paper I) in a kinetic 
approach (e.g. [8]). 

The kinetic model prediction for the peak value of the expected 7-ray flux is 
not very different from that obtained in the simplified models [5] , even though 
this difference is not unimportant. The main reason is that this peak value is 
mainly determined by the fraction of the explosion energy that is converted 
into CR energy, i.e. by the efficiency of CR acceleration, which is not strongly 
dependent on the model used. 

There are more important differences in the time variation of the predicted 
7-ray fluxes. Kinetic theory revealed a much more effective CR and therefore 
7-ray production during the free expansion phase of the SNR, and a more rapid 
decrease of the 7-ray flux after reaching its peak value during the subsequent 
Sedov phase, due to the effect of different spatial distributions of the gas and 
the CRs inside the SNR (the so-called overlapping effect) that had not been 
taken into account in [5]. 

In Paper I we have considered only one particular set of possible physical 
parameters. In the present calculations we demonstrate how the expected TeV- 
energy 7-ray flux from SNRs depends upon the ISM density, the ejected mass 
and the CR diffusion coefficient. 



2 Results and Discussion 

We use here the full time- dependent, kinetic model for particle acceleration 
in SNRs which selfconsistently describes diffusive shock acceleration of CRs, 
taking account the nonlinear CR backreaction on the structure and dynamical 
evolution of the expanding spherical supernova (SN) shock [8]. CRs naturally 
originate from a suprathermal gas particle population at the shock front. This 
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means that after shock heating the most energetic gas particles become in- 
volved, i.e. "injected", into the acceleration process. 

Detailed investigations of CR acceleration and SNR evolution in the uniform 
ISM have demonstrated important features of this process [8]. For a wide 
range of possible injection rates the CR acceleration efficiency is very high and 
almost independent of the injection rate. This is why we use here the particular 
value of the injection parameter r\ = 10~ 4 , which denotes the fraction of gas 
particles involved in the acceleration. This value of i] provides an injection 
rate which is more than an order of magnitude lower compared with results 
of collisionless shock plasma simulations [9,10] and analytical injection theory 
[11], and with an injection rate that corresponds to the kinetic Monte Carlo- 
model (e.g. [12]) for purely parallel shocks. Our chosen value of 77 effectively 
takes the influence of the shock obliquity into account: according to [12,13], 
already at angles 9 45° between the upstream magnetic field and the shock 
normal the injection rate is about an order of magnitude smaller than in the 
purely parallel shock case. 

As such the process of CR acceleration and associated 7-ray production in 
SNRs is detailed in Paper I. Here we present only calculations of the expected 
TeV 7-ray emission, measureable by the imaging Cherenkov technique. In or- 
der to illustrate the sensitivity of the 7-ray production to relevant physical 
parameters of the ISM, we present in Fig. 1 calculations of the expected in- 
tegral 7-ray flux F 7 (> 1 TeV), normalized to the distance d — 1 kpc, for 
eight different parameter sets. Four of them correspond to physical parame- 
ters typical for the case of a SN la: SN explosion energy E sn = 10 51 erg, ejecta 
mass M e j = 1.4M Q , and the value k = 7 of the parameter k which describes 
the ejecta velocity distribution; four other parameters correspond to the core 
collapse SN lb /SN II cases: E sn = 10 51 erg, M ej = 10 M , k = 10. 

Note that an essential part of the volume of our Galaxy is occupied by a 
rarefied, so-called hot ISM phase, with a density that is about two orders of 
magnitude lower than another representative phase, the so-called warm ISM 
phase, with hydrogen number density Nh = 0.3 cm -3 . The kinetic theory 
results imply that CRs are effectively produced also by SNRs in the hot ISM. 
But even simple estimates show that in this case the expected 7r°-decay 7-ray 
flux is far below practical detection possibilities, due to the extremely low 
ISM density (there might be detectable Inverse Compton emission, which is 
however not considered here). Therefore, in both cases above, the calculations 
were performed for an ISM hydrogen number density N H = 0.3 cm" 3 and 30 
cm -3 . The latter case can model a SNR evolving inside a rather dense cloud. 

SN 1006, the remnant of a SN la, is the only SNR where TeV 7-ray emis- 
sion was definitely claimed to have been detected up to now [14]. It exhibits 
also strong X-ray synchrotron emission with a spectrum that implies a maxi- 
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Fig. 1. Integral TeV 7-ray flux, normalized to the distance d = 1 kpc, as a function 
of time. Thick lines correspond to the Bohm diffusion coefficient kb(p), thin lines 
correspond to k = IOkb- Four different cases correspond to the following 3-param- 
eter set: M ej /M & , N H /(1 cm" 3 ), t /(l yr)): 1.4, 0.3, 368 (full lines); 1.4, 30, 79 
(dashed lines); 10, 0.3, 1897 (dash-dotted lines); 10, 30, 409 (dotted lines). 

mum electron energy e max > 10 TeV. According to the calculations in [15], this 
should also be the maximum energy of the proton component which implies 
that CR diffusion is about ten times less effective than Bohm diffusion. There- 
fore we shall also consider two different values of the CR diffusion coefficient. 
The first corresponds to the ordinary Bohm-type diffusion coefficient k = Kb 
(e.g. Paper I); in the second case we use a ten times higher value, k = IOkb- 
Note that different values of k influence only the value of the maximum CR 
energy e max (or maximum momentum p max = tmax/c, where c is the speed of 
light), reached during SNR evolution. In all cases we use the same standard 
value of the ISM magnetic field, B = 5 fiG. 



To the extent that the CR energy e ~ 10e 7 ~ 10 TeV , responsible for 7-rays 
with energy e 7 = 1 TeV, is close to e max , k essentially influences the expected 
7-ray flux F 7 (> 1 TeV): in the case e max < 10 TeV F 7 (> 1 TeV) should 
be essentially lower than at e max ^10 TeV. Therefore we present in Fig. 2 the 
calculated value of the maximum CR momentum p max as a function of time. By 
definition p max is the CR momentum, where the overall momentum spectrum 
N(p, t) of CRs, accelerated up to time t, deviates from the dependence p~ 2 by 
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Fig. 2. The maximum CR momentum as a function of time for the same cases as 
in Fig. 1. 

a factor of e: 



For the sake of convenience, in Figs. 1 and 2 time is measured in units of to, 
where t is related to the sweep-up radius R and the mean initial ejecta speed 
V as follows: 



Here po — 1.4 7V#m is the density of the ISM which contains 10% helium, and 
m is the proton mass. 

The maximum momentum p max is determined by geometrical factors [16]. 
The CR spectrum starts to deviate significantly from its power law form at 
momenta p ~ p ma x, where the shock becomes too slow - and too small in 
radial extent - to fill the upstream region with CRs with a number density 
that is sufficient to provide a power law spectrum. The shock size R s and its 
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speed V s are the most relevant parameters which influence the value of p m ax- 
In the free expansion phase, t < to, p max R S V S is an increasing function 
of time [16]. During the Sedov phase, t^t , p ma x remains almost constant, 
because the shock weakens at this stage, the value of R S V S oc t -1 / 5 decreases 
with time, and the shock produces CR particles only up to the momentum 
p m oc R S V S which is a decreasing function of time. During this phase the value 
of Pmax is determined by the CRs accelerated at the end of the free expansion 
phase. For t ^> t the influence of the shock on these highest energy CRs 
(the so-called escaping particles) becomes almost negligible, and their energy 
remains nearly constant. The calculated results correspond very well to the 
expected dependence [16]: 

p max oc eH 2 M- 1/6 N h 1/3 /(k b /k). (3) 



However this relation was derived in the test particle approximation and does 
not take into account the nonlinear effects which also influence p m ax- As one 
can see from Fig. 2, the exact value of p ma x in the Sedov phase is slightly 
higher for a larger ejected mass M e j, contrary to relation (3): the modified 
SN shock produces high energy CRs more efficiently and this increases their 
maximum momentum by a factor of a few compared with the case of an 
unmodified shock. For example, in the case k = Kb, Nh = 0.3 cm -3 , and 
M e j = 10M Q , the shock becomes significantly modified already at the end of 
the free expansion phase. The shock compression ratio is a = 5.5 at t — to , 
and reaches its maximum value 5.8 at t — 2t . Therefore the maximum CR 
momentum reaches p max w 2 x 10 5 mc soon after the beginning of the Sedov 
phase (see Fig. 2). It is slightly larger than the estimate 1.6 x 10 5 mc predicted 
by test particle theory [16]. Due to the higher mean ejecta speed Vo in the case 
M e j = IAMq, the shock reaches peak modification significantly later: it is only 
slightly modified (a = 4.5) at t — to and the maximum shock compression ratio 
a = 5.8 is reached at t — 10t . At this stage, when the shock produces CRs 
most effectively, the product R S V S has a lower value than at t — to. Therefore 
the maximum value p ma x — 1-8 x 10 5 mc reached at t — 10t is slightly smaller 
than the test particle prediction 2.2 x 10 5 mc. The effect described gives in 
the Sedov phase a maximum CR momentum which is higher for larger ejecta 
mass, opposite to the test particle prediction (3), although the exact value of 
Pmax does not deviate essentially from the test particle prediction. 

As can be seen from Fig. 2, for Nh = 0.3 cm -3 the maximum CR momentum 
remains above the critical level p m axif) ~ 10 4 mc at least for t/t ~ 0.1 even for a 
large diffusion coefficient k = 10 k,b- Therefore the value of k does not influence 
so much the expected flux F 7 (> 1 TeV). In this case the predicted flux F 1 is 
roughly proportional to the ISM density over the whole SNR evolution. The 
time profile of this integral 7-ray flux above 1 TeV is denoted by F 7 (t/t ) for 
short, and is roughly the same in all four cases. Note that the essential part 
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of the dependence of F 7 (t/to) on M e j and Nh is already contained in the time 
unit to. The remaining part of the dependence on M e j is the dependence of 
the time t m (at which the peak value of the 7-ray flux is reached) upon M e j. 
As one can see from Fig. 1 (again using N H = 0.3 cm -3 ), in the case k = k b 
t m /t ps 2.5 for M ej = 1OM , and t m /t ^ 10 for M ej = 1.4M . In the case 
k = IOkb, t m /t ~ 6.5 for M ei = 1OM , and * m /t ~ 32 for M ej = 1.4M . 
The reason is that one of the most relevant factors determining the value of the 
7-ray flux is the CR energy content in the SNR. It reaches its peak value in the 
Sedov phase, when the shock speed V s drops to some critical value V m , which 
is a function of the ISM parameters (see [7] for details). Since t oc M^j 6 , we 

have the relation t m /t oc M~- 5 ^ 6 . This corresponds satisfactorily to the results 
represented in Fig. I. 

Only in the case of a dense ISM, N H = 30 cm" 3 , does the increase of CR 
diffusion coefficient lead to a drop of the CR maximum momentum below the 
critical level p max < 10 4 mc, that causes F 7 (> 1 TeV) to drop significantly 
compared with its value in the case k = k,b- In this case TeV 7-rays are 
produced by CRs which belong to the exponential part of their spectrum 
P^Pmax- Therefore the 7-ray emissivity drops exponentially with increasing 
energy e 7 . Note that this significant influence of the CR diffusion coefficient 
on the expected 7-ray flux takes place only in the energy range e 7 > 1 TeV for 
this density. For considerably lower energies e 7 < 1 TeV the 7-ray production 
is due to correspondingly lower-energy particles which belong to the power 
law part of CR spectrum, and thus the spectrum remains insensitive to the 
value of the CR diffusion coefficient. As can be seen from Fig. 1, the expected 
TeV 7-ray flux remains almost unchanged despite of an increase of the ISM 
density from Nh = 0.3 cm" 3 up to 30 cm" 3 in the case k = 10/c.b. This implies 
that the effect of the /t-increase is about 100. 

The highest energy CRs with momenta p^p ma x become almost insensitive 
to the shock influence at the late Sedov phase t ^> to [8]. They fill a volume 
almost uniformly whose size increases with time at this stage according to the 
diffusive law R oc y/Kt. Hence we have the situation when CRs with a number 
density n(e) oc i?" 3 interact with a progressively increasing amount of gas 
M oc i? 3 , resulting in an almost constant TeV 7-ray flux F 7 oc Mn (see the 
curves, which correspond to k — 10k b and N H = 30 cm" 3 , in Fig. 1). 

Our calculations show that the efficiency of TeV 7-ray production by shock ac- 
celerated CRs in a SNR is characterized by a critical value of the CR diffusion 
coefficient, that can be represented in the form 



where the value K ~ 10 is expected to be only slightly dependent on the 
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injection rate. For CR diffusion coefficients significantly larger than n C rit, we 
expect maximum CR energies e max < 10 TeV, which would lead to a sizeable 
decrease of the TeV 7-ray flux below a detectable level at all ISM densities 
considered. 

The critical value of the CR diffusion coefficient exists for an arbitrary 7-ray 
energy e 7 . It is clear from the above that it is determined by the expression 
(4) with 



This simply shows that the lower 7-ray energy e 7 is, the wider is the range 
of the CR diffusion coefficient which allows efficient production of 7-rays with 
energy e 7 . 



3 Summary 

We have studied here how the background ISM density and the CR diffusion 
coefficient influence the expected 7r°-decay 7-ray production in SNRs in the 
TeV range. Our calculations show that a peak TeV 7-ray flux F m , normalized 
to a distance of d — 1 kpc, of about 10~ 10 photons cm" 2 s" 1 is reached at 
t m ~ 3to, for an ISM number density Nh = 0.3 cm -3 . In the Bohm limit 
for CR diffusion the ratio F m /Nn is independent of ISM density in the range 
Nh = 0.3-f-30 cm" 3 . If the 7-ray flux is written in the form _F 7 (t/t ), then the 
ejected mass M e j influences only the peak time t m , with t m /t oc M e j . 

If CR scattering is ten times less efficient near the shock front, then the max- 
imum CR energy decreases below the critical value 10 TeV which makes the 
7-ray production almost independent of the ISM density, at a level that cor- 
responds to N H = 0.3 cm -3 in the Bohm limiting case. 

If the CR diffusion coefficient is significantly larger than K crit (see eq.(4)), 
then the expected TeV 7-ray flux will drop considerably below the presently 
detectable level of 10~ 12 photons cm" 2 s" 1 . Note that in this case SNRs would 
also be hardly considered as the sources of the Galactic CRs. 
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